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Baryonic Operators for Lattice Simulations 

LHP Collaboration: R. Edwards^, R. Fiebig'^, G. Fleming^, U.M. Heller^ C. Morningstar'^ *, 
D. Richards^, I. Sato°, S. Wallace^ 

^Thomas Jefferson National Accelerator Facility, Newport News, VA 23606, USA 
^'Physics Department, Florida International University, Miami, FL 33199, USA 
American Physical Society, One Research Road, Ridge, NY 11961-9000, USA 
'^Department of Physics, Carnegie Mellon University, Pittsburgh, PA 15213, USA 
^Department of Physics, University of Maryland, College Park, MD 20742, USA 

The construction of baryonic operators for determining the A^* excitation spectrum is discussed. The operators 
are designed with one eye towards maximizing overlaps with the low-lying states of interest, and the other eye 
towards minimizing the number of sources needed in computing the required quark propagators. Issues related 
to spin identification are outlined. Although we focus on tri-quark baryon operators, the construction method is 
applicable to both mesons and penta-quark operators. 



1. OVERVIEW 

One major goal of the Lattice Hadron Physics 
Collaboration (LHPC) is to calculate a significant 
portion of the low-lying spectrum of baryon reso- 
nances from Monte Carlo estimates of Euclidean- 
space path integrals in lattice QCD. To accom- 
plish this goal, large sets of operators which can 
create the baryon states of interest at a given time 
t must be designed. Our approach to constructing 
such operators is outlined in this poster. 

Baryon states are identified by their momen- 
tum p, intrinsic (half-integral) spin J, projection 
A of this spin onto some axis, parity P = ±1, and 
quark flavor content (isospin, strangeness, etc.). 
Since we are interested first in the masses of these 
states, we restrict our attention to the p = 
sector, so our operators must be invariant under 
all spatial translations allowed on a cubic lattice. 
The little group of all symmetry transformations 
which leave p — invariant is the octahedral 
point group Oh, so our operators may be classi- 
fied using the spinorial irreducible representations 
(irreps) of Oh ■ There are four two-dimensional ir- 
reps Gig, Giu, G2g, G2u and two four-dimensional 



representations Hg and H„ |1I2| . The continuum- 
limit spins J of our states must be deduced by 
examining degeneracy patterns across the differ- 
ent Oh irreps. 

We use a two-step approach to construct the 
operators. First, we design elemental operators 
[t, x) having the appropriate flavor structure 
characterized by isospin, strangeness, etc., and 
color structure constrained by gauge invariance. 
In the second step, we apply group-theoretical 
projections to obtain operators which transform 
irreducibly under all lattice rotation and reflec- 
tion symmetries: 

Bt^'^it, X) = A^i^(A)*(i?)[/^ijf (t, ^)ul (1) 

where A refers to an Oh irrep, A is the irrep row, 
is the number of elements in Oh, is the 
dimension of the A irrep, DmXi^R) is a A repre- 
sentation matrix corresponding to group element 
R, and Ur is the quantum operator which imple- 
ments the symmetry operations. We also impose 
translation invariance for zero momentum states: 
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Note that in order to capture the spectrum^ we 
plan to calculate correlation matrices 



C^^^it) = (0| TBt^^m'^^^iQ) |0). (3) 



Matrices of correlation functions facilitate the re- 
liable extraction of the masses of excited states 
since variational techniques, possibly combined 
with Bayesian methods, can be exploited. The 
operators are designed with one eye towards max- 
imizing overlaps with the low-lying states of in- 
terest, and the other eye towards minimizing the 
number of sources needed to compute the re- 
quired quark propagators. Throughout this work, 
we work in the approximation to„ = m^, where 
niu and nid are the masses of the u and d quarks, 
respectively, and we ignore electromagnetic ef- 
fects. Hence, isospin is a good quantum number. 

2. DESIGN DETAILS 

Gauge invariance and isospin dictate the color 
and flavor structure of the elemental operators in 
the first stage of the operator construction. Since 
we intend to include all allowed operators within a 
given isospin channel, we do not need to consider 
explicit S'?7(3)-flavor combinations. As long as 
we are careful to respect isospin, the variational 
method takes care of finding the preferred flavor 
combinations for us. 

Both link variable and quark field smear- 
ings are important for reducing the couplings of 
our operators to higher-lying states. Familiar 
APE smearing^ of the spatial links can be em- 
ployed, and the three-dimensional covariant spa- 
tial Laplacian is used to smear the quark fields: 



i;{x) = (l + eA)"« 7p{x), 



(4) 



where g and integer Ug are tunable parameters 
and the three-dimensional gauge-covariant Lapla- 
cian is defined by 

K^{x)^ (c/fc(a;)V(a; + fc)-^(x)), (5) 

fc=±l,±2,±3 

with U^{x) denoting a smeared link variable. 
Note that the square of the smeared field is equal 
to zero, similar to the simple Grassmann field. 



Our three-quark elemental operators are con- 
structed using the following building blocks: 

xS,„(x) ^ (AV(x))^^^, (6) 

eZi(-) ^ {D^K-i^ix)) (7) 

where A is a flavor index, a is a color index, a is a 
Dirac spin index, and the p-link gauge-covariant 
displacement operator is defined by 

b\^^0{x) = U,ix) . . . U,{x+ip-l)j)0{x+pj).i8) 

Our three-quark elemental operators have the fol- 
lowing forms: 



^ABC ^abc XAaa XBb0 Xcc-f^ 



rii 712 f^^i 
■abc XAaa XBb0 ^C. 



-'ABC ^abc XAaa ^BbB 



Bbl3 ^Cc-f ' 



(9) 
(10) 

(11) 



where ni, n2, n3,p,pi,p2 are positive integers, 
j,k = ±1,±2,±3 are spatial directions, a,/3, 7 
are Dirac spin indices, A, _B, C are flavors, and 
a, 5, c are colors. The Levi-Civita symbol £abc 
specifles the color structure and the 4>^bc specify 
the flavor combinations. Since baryon resonances 
are expected to be large objects, extended oper- 
ators are crucial. Note that the different powers 
of the spatial Laplacian A are used to build up 
radial structure, while the displacement opera- 
tor d'"^^ is used to incorporate orbital structure. 
A large number of baryons can be studied using 
a somewhat small number of quark propagator 
sources: ni^{na -I- 3npnb), where Up is the num- 
ber of displacement lengths p, Ua is the number 
of powers of A, rii, is the number of powers of A 
combined with a displacement, and is number 
of quark masses used. 

In the second stage of the operator construc- 
tion, the group-theoretical projections in Eq. 
are applied. Eq. (^) requires explicit representa- 
tion matrices for every group element. Represen- 
tation matrices for all allowed proper rotations 
can be generated from the matrices for C^y and 
C4Z, the rotations by 7r/2 about the y- and z- 
axes, respectively. The explicit matrices we use 
for these generators are given by 
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Table 1 

Continuum limit spin identification: occurrences 
rip of O irrep T in subduction of SU{2) irrep J. 
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1 

71 



l-i 
l + i 
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-V3 


x/3 


-1 ' 
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x/3 


-1 


-1 
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V3 
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-V3 
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V3 


n/3 


1 . 
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Inclusion of spatial inversion is straightforward. 

The projections in Eq. are carried out us- 
ing Maple. The resulting operators transform ir- 
reducibly under Oh- In the continuum limit, the 
spin J is identified using Tabled For example, to 
identify a spin-^ baryon, it must occur in the Gi 
irrep and accompanying degenerate levels in the 
G2 and H irreps must not occur. A level occur- 
ring at the same a —> Q mass in all three Gi, G2, 
and H irreps can be identified as a spin-| baryon. 

To illustrate our operator construction, con- 
sider an explicit example. Let 



af3-y 



(12) 



ignoring smearing and powers of the Laplacian 
to simplify matters and using u, d to indicate fla- 
vor. The above operators have isospin / = 1/2, 
isospin projection 1^ = 1/2, and strangeness 
5 = 0. There are twenty independent opera- 
tors due to the constraints ^ap-^ + ^-yiia = and 
+ ^/37a + ^laf} = 0. Linear combinations 
which transform irreducibly under Oh are listed 
in Table El 



Table 2 

Combinations of the operators ^ap-^ in Eq. H12|l 
which transform irreducibly under Oh for the 
Dirac-Pauli representation of the 7-matrices. 



Irrep 


Row 


Operators 


Gig 


1 


*221 


Gig 


2 


*112 


Gig 


1 


$324 


Gig 


2 


-*314 


Gig 


1 


2$44i + 2$234 - *324 


Gig 


2 


2$332 + 2$i34 - $314 


Gi„ 


1 


*443 


Gijj 


2 


$334 


Gijj 


1 


$124 — *i>214 


Giu 


2 


— *I'l23 + *I'213 


Gi„ 


1 


2$223 — *i'l24 — $214 


Gi„ 


2 


2<I>114 - $123 - $213 


Ha 


1 


V3 $442 


Hg 


2 


-$441 + 2$234 - $324 


Ha 


3 


$332 - 2$i34 + $314 


Hg 


4 


-V3 $331 


Hu 


1 


V3 $224 


Hu 


2 


— $223 — $124 — $214 


Hu 


3 


$114 + $123 + $213 


Hu 


4 


-\/3 $113 



Although we focused on tri-quark baryon oper- 
ators, our approach to constructing hadronic op- 
erators is applicable to mesons and penta-quark 
operators. By including E and B fields into our 
operators, hadrons bound by an excited gluon 
field can also be studied. This work was sup- 
ported by the U.S. National Science Foundation 
under Awards PHY-0099450 and PHY-0300065, 
and by the U.S. Department of Energy under 
contracts DE-AC05-84ER40150 and DE-FG02- 
93ER-40762. 
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